function main() {
	let scale=255.0/4000 ;
	let ds=pe.Dataset('modis/mod09cmg', 20230901000000);
	let ds1 = ds.forEachPixel( function(v,i) {
		r=v[0];//modis band1 for red
		g=v[3];//modis band4 for green
		b=v[2];//modis band3 for blue
		if( r>0 ) {
			r=r*scale;
			g=g*scale;
			b=b*scale;
		}else{
			r=g=b=0;
		}
		return [r,g,b] ;
	}) ;
	let geomPoint=pe.geomByWkt("POINT (123.046875 36.738884)").buffer(2.0);
	let geomMPoints=pe.geomByWkt("GEOMETRYCOLLECTION(POINT (87.1875 48.922499),"
		+" POINT (106.171875 42.55308), POINT (120.410156 53.120405))").buffer(1);
	let geomLine=pe.geomByWkt("LINESTRING (-20.3764 -18.3506, -75.36948 -18.34094,"
		+" -80.36277 -18.33377, -85.35619 -18.32807, -90.34925 -18.32336, -95.34269"
		+" -18.31949, -100.33627 -18.31624, -105.32741 -18.3125, -110.32558 -18.31201,"
		+" -116.3764 -18.3506 )").buffer(1);
	let geomPoly=pe.geomByWkt("POLYGON ((-75.234375 10.487812, -82.96875 -5.441022,"
		+" -72.070313 -18.479609, -53.085938 -22.755921, -40.605469 -18.979026,"
		+" -37.617188 -12.726084, -33.925781 -5.965754, -40.078125 -1.933227,"
		+" -50.097656 1.230374, -53.789063 7.885147,"
		+" -61.171875 10.487812, -71.894531 12.21118, -75.234375 10.487812))").buffer(1);
	ds1=ds1.burn(geomPoint,22) ;
	ds1=ds1.burn(geomMPoints,22) ;
	ds1=ds1.burn(geomLine,22) ;
	ds1=ds1.burn(geomPoly,22) ;
	return ds1;
}